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Abstract — Low-density parity check (LDPC) codes are an 
important class of codes with many applications. Two algebraic 
methods for constructing regular LDPC codes are derived - one 
based on nonprimitive narrow-sense BCH codes and the other 
directly based on cyclotomic cosets. The constructed codes have 
high rates and are free of cycles of length four; consequently, they 
can be decoded using standard iterative decoding algorithms. 
The exact dimension and bounds for the minimum distance and 
stopping distance are derived. These constructed codes can be 
used to derive quantum error-correcting codes. 

Index Terms — LDPC Codes, BCH Codes, Channel Coding, 
Performance and iterative decoding, quantum BCH codes. 

I. Introduction 

Bose-Chaudhuri-Hochquenghem (BCH) codes are an inter- 
esting class of linear codes that has been investigated for 
nearly half of a century. These types of codes have a rich 
algebraic structure. BCH codes with parameters [n, k,d> S] q 
are interesting because one can choose their dimension k 
and minimum distance d once given their design distance 6 
and length n over a finite field with q elements. A linear 
code defined by a generator polynomial g(x) has dimension 
k = n — deg(g(x)) and rate k/n. It is not easy to show the 
dimension of nonprimitive BCH codes over higher finite fields. 
In [3], [4], we have given an explicit formula for the dimension 
of these codes if their deigned distance S is less than a constant 

^max- 

Low-density parity check (LDPC) codes are a capacity- 
approaching class of codes that were first described in a 
seminal work by Gallager [9]. Tanner in [20] rediscovered 
LDPC codes using a graphical interpretation. A regular (p, A) 
LDPC code is measured by the weights of its columns p and 
rows A. Iterative decoding algorithms of LDPC and turbo 
codes highlighted the importance of these classes of codes 
for communication and storage channels. Furthermore, these 
codes are practical and have been used in many beneficial 
applications [6], [12]. In contrast to BCH and Reed-Solomon 
(RS) cyclic codes, LDPC cyclic codes with sparse parity check 
matrices are customarily constructed by a computer search. In 
practice, LDPC codes can achieve higher performance and 
better error correction capabilities than many other codes, 
because they have efficient iterative decoding algorithms, such 
as the product-sum algorithm [12]— [14], [21]. Some BCH 
codes turned out to be LDPC cyclic codes as well; for example, 



a [15, 7] BCH code is also an LDPC code with a minimum 
distance five. 

Regular and irregular LDPC codes have been constructed 
based on algebraic and random approaches [12], [18], and 
references therein. Liva et al. [13] presented a survey of the 
previous work done on algebraic constructions of LDPC codes 
based on finite geometry, elements of finite fields, and RS 
codes. Yi et al. [22] gave a construction for LDPC codes, 
based on binary narrow-sense primitive BCH codes, and their 
method is free of cycles of length four. Furthermore, a good 
construction of LDPC codes should have a girth of the Tanner 
graph, of at least 6 [12], [13]. One might wonder how do 
the rates and minimum distance of BCH codes compare to 
LDPC codes? Do self-orthogonal BCH codes give raise to self- 
orthogonal LDPC codes as well under the condition 8 < S max . 
We show that how to derive LDPC codes from nonprimitive 
BCH codes. 

One way to measure the decoding performance of linear 
codes is by computing their minimum distance d m i n . The 
performance of low-density parity check codes under iterative 
decoding can also be gauged by measuring their stopping sets 
S and stopping distance s, which is the size of the smallest 
stopping set [16], [17]. For any given parity check matrix H 
of an LDPC code C, one can obtain the Tanner graph G of this 
code and computes the stopping sets. Hence, s is a property 
of H, while d m i n is a property of C. The minimum distance is 
also bounded by d m i n > s. BCH codes are decoded invertible 
matrices such as Berkcampe messay method, LDPC codes ar 
decoded using iterative decoding and Belief propagation (BP) 
algorithms. 

In this paper, we give a series of regular LDPC and Quasi- 
cyclic (QC)-LDPC code constructions based on non-primitive 
narrow-sense BCH codes and elements of cyclotomic cosets. 
The constructions are called Type-I and Type-II regular 
LDPC codes. The algebraic structures of these codes help us 
to predict additional properties of these codes. Hence, The 
constructed codes have the following characteristics: 

i) Two classes of regular LDPC codes are constructed that 
have high rates and free of cycles of length four. Their 
properties can be analyzed easily. 

ii) The exact dimension is computed and the minimum 
distance is bounded for the constructed codes. Also, the 
stopping sets and stopping distance can be determined 
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from the structure of the parity check matrices. 
The motivation for our work is to construct Algebraic 
regular LDPC codes that can be used to derive quantum 
error-correcting codes. Alternatively, they can also be used for 
wireless communication channels. Someone will argue about 
the performance and usefulness of the constructed regular 
LDPC codes in comparison to irregular LDPC codes. Our 
first motivation is to derive quantum LDPC codes based on 
nonprimitive BCH codes. Hence, the constructed LDPC-BCH 
codes can be used to derive classes of symmetric quantum 
codes [l]-[3], [5], [15] and asymmetric quantum codes [8], 
[19]. The literature lacks many constructions of algebraic 
quantum LDPC codes, see for example [1], [15] and references 
therein. 

II. Constructing LDPC Codes 

Let F g denote a finite field of characteristic p with q 
elements. Recall that the set F* = F g \ {0} of nonzero field 
elements is a multiplicative cyclic group of order q — 1. A 
generator of this cyclic group is called a primitive element of 
the finite field F„. 



A. Definitions 

Let n be a positive integer such that gcd(n, q) = 1 and 
/i = q m — 1, where m = ord„(g) is the 



Q 



|m/2j 



< n < 



multipicative order of q modulo n. 

Let a denote a fixed primitive element of F g m. Define a 
map z from F* m to F% such that all entries of z(a z ) are 
equal to except at position i, where it is equal to 1. For 
example, z(a 2 ) = (0,1,0,..., 0). We call z(a k ) the location 
(or characteristic) vector of or. We can define the location 
vector z(a J+J+1 ) as the right cyclic shift of the location vector 
z(a l+ - J ), for < j < fi — 1, and the power is taken module 
fi. 

Definition 1: We can define a map A that associates to an 



element F* 



a circulant matrix in F% Xfi 



by 



z(a') 
z(a i+1 ) 

\ z(a l +^ 1 ) ) 



(1) 



By construction, A(a k ) contains a 1 in every row and column. 

For instance, A(or ) is the identity matrix of size fix fi, and 
A(a 2 ) is the shift matrix 



A(a 2 ) 



I 1 
1 



V 1 







(2) 



We will use the map A to associate to a parity check matrix 
H = (hij) in (F* m ) axb the (larger and binary) parity check 



matrix H = (A(hij)) in F 



The matrices A(hij)'s are 



fix fi circulant permutation matrices based on some primitive 
elements hij as shown in Definition Q] 



B. Regular LDPC Codes 

A low-density parity check code (or LPDC short) is a binary 
block code that has a parity check matrix H in which each row 
(and each column) is sparse. An LDPC code is called regular 
with parameters (p, A) if it has a sparse parity check matrix 
H in which each row has p nonzero entries and each column 
has A nonzero entries. 

A regular LDPC code defined by a parity check matrix H is 
said to satisfy the row-column condition if and only if any two 
rows (or, equivalently, any two columns) of H have at most 
one position of a nonzero entry in common. The row-column 
condition ensures that the Tanner graph does not have cycles 
of length four. 

A Tanner graph of a binary code with a parity check matrix 
H = (hij) is a graph with vertex set V U C that has one 
vertex in V for each column of H and one vertex in C for 
each row in H, and there is an edge between two vertices i and 
j if and only if hij ^ 0. Thus, the Tanner graph is a bipartite 
graph. The vertices in V are called the variable nodes, and 
the vertices in C are called the check nodes. We refer to d(vi) 
and d(cj) as the degrees of variable node v,i and check node 
Cj respectively. 

Two values used to measure the performance of the decod- 
ing algorithms of LDPC codes are: girth of a Tanner graph 
and stopping sets. The minimum stopping set is analogous to 
the minimum Hamming distance of linear block codes. 

Definition 2 (Grith of a Tanner graph): The girth g of the 
Tanner graph is the length of its shortest cycle (minimum 
cycle). 

A Tanner graph with large girth is desirable, as iterative 
decoding converges faster for graphs with large girth. 

Definition 3 (Stopping set): A stopping set S of a Tanner 
graph is a subset of the variable nodes V such that each vertex 
in the neighbors of S is connected to at least two nodes in S. 

The stopping distance is the size of the smallest stopping 
set. The stopping distance determines the number of cor- 
rectable erasures by an iterative decoding algorithm, see [7], 
[16], [17]. 

Definition 4 (Stopping distance): The stopping distance of 
the parity check matrix H can be defined as the largest integer 
s(H) such that every set of at most (s(H) — 1) columns of H 
contains at least one row of weight one, see [17]. 
The stopping ratio a of the Tanner graph of a code of length 
n is defined by s over the code length. 

The minimum Hamming distance is a property of the code 
used to measure its performance for maximum-likelihood 
decoding, while the stopping distance is a property of the 
parity check matrix H or the Tanner graph G of a specific 
code. Hence, it varies for different choices of H for the same 
code C. The stopping distance s(H) gives a lower bound of 
the minimum distance of the code C defined by H, namely 



s(H) < d n 



(3) 



It has been shown that finding the stopping sets of minimum 
cardinality is an NP-hard problem, since the minimum-set 
vertex covering problem can be reduced to it [11]. 
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III. LDPC Codes based on BCH Codes 

In this section we give two constructions of LDPC codes 
derived from nonprimitive BCH codes, and from elements 
of cyclotomic cosets. In [22], the authors derived a class of 
regular LDPC codes from primitive BCH codes but they did 
not prove that the construction has free of cycles of length four 
in the Tanner graph. In fact, we will show that not all primitive 
BCH codes can be used to construct LDPC with cycles greater 
than or equal to six in their Tanner graphs. Our construction 
is free of cycles of length four if the BCH codes are chosen 
with prime lengthes as proved in Lemma UJ in addition the 
stopping distance is computed. Furthermore, We are able to 
derive a formula for the dimension of the constructed LDPC 
codes as given in Theorem [9] We also infer the dimension 
and cyclotomic coset structure of the BCH codes based on 
our previous results in [3], [4]. 

We keep the definitions of the previous section. Let q be a 
power of a prime and n a positive integer such that gcd(g, n) = 
1. Recall that the cyclotomic coset C x modulo n is defined as 



C x = {xq l mod n \ i <G Z, i > 0}. 



(4) 



Let m be the multiplicative order of q modulo n. Let a 
be a primitive element in F g ™ . A nonprimitive narrow-sense 
BCH code C of designed distance 8 and length n over F q is a 
cyclic code with a generator monic polynomial g(x) that has 



a, a 2 , . . . , a 5 1 as zeros, 



5-1 



(5) 



Thus, c is a codeword in C if and only if c(a) = c(a 2 ) = 
. . . = c(a <5_1 ) = 0. The parity check matrix of this code can 
be defined as 



H, 



belt 



a 



a 5 - 1 a 2 ^ 



a 

T,2(n-1) 



(5-l)(n-l) 



(6) 



We note the following fact about the cardinality of cyclo- 
tomic cosets. 

Lemma 5: Let n be a positive integer and q be a power of 
a prime, such that gcd(n, q) — 1 and gL m / 2 J < n < q™ — 1, 
where m = ord n (q). The cyclotomic coset C x = {xqi mod 
n | < j < m} has a cardinality of m for all x in the range 
1 < x < nq^ m ^/{q m - 1). 

Proof: See [3, Lemma 8]. ■ 

Therefore, all cyclotomic cosets have the same size m if 
their range is bounded by a certain value. This lemma enables 
one to determine the dimension in closed form for BCH code 
of small designed distance [3], [4]. In fact, we show the 
dimension of nonprimitve BCH codes over F q . 

Theorem 6: Let q be a prime power and gcd(n, q) = 1, 
with ord n (q) = m. Then a narrow-sense BCH code of length 



Jm/2J 



the range 2 < 8 < <5 max = min{ [nq^ m ^/(q m - l)J,n}, has 
dimension of 



-m|-(*-l)(l-l/g)l. 



(7) 



Proof: See [3, Theorem 10]. ■ 
Based on these two observations, we can construct regular 
LDPC codes from BCH codes with a known dimension and 
cyclotomic coset size. 

A. Type-I Construction 

In this construction, we use the parity check matrix of a 
nonprimitive narrow-sense BCH code over F 9 to define the 
parity check matrix of a regular LDPC over F2. 

Consider the narrow-sense BCH code of prime length 
q[m/2\ <; n < qm _ ^ over jf^ w j t j 1 designed distance 5 

and ord n (q) = m. We use the fact that there must be some 
primes in the integer range (qL m / 2 J g m _ 1) in f ac t ; there 
must exist a prime between x and 2x for some integer x, 
in which it ensures existence primes in the given interval. A 
parity check matrix H of an LDPC code can be obtained by 
applying the map A in Equation (Q~|) to each entry of the parity 
check matrix (O of this BCH code, 



H 



A(l) 
A(l) 



A(a) 
A(a 2 ) 



A(a 2 ) 
A(a 4 ) 



A(a™- 1 ) 
A(a 2 ("" 1 )) 



(8) 



A(l) A(a 5 - 1 ) A{a 2 ^-^) 



A(a 



(S-l)(n-iy 



The matrix H is of size (i-lj/ixw/j and by construction it 
has the following properties: 

• Every column has a weight of 8 — 1. 

• Every row has a weight of n. 

The matrix H of size (5—l)fixnfi has a weight of p = 5—1 
in every column, and a weight of A = n in every row. The 
null space of the matrix H defines a (p, A) LDPC code with 
a high rate for a small designed distance 8 as we will show. 
The minimum distance of the BCH code is bounded by 



> 



1, odd 8; 

2, even 8. 



(9) 



Also, the minimum distance of the LDPC codes is bounded by 
drain- Now, we will show that in general regular (p, A) LDPC 
codes derived from primitive BCH codes of length n are not 
free of cycles of length four as claimed in [22]. 

Lemma 7: The Tanner graph of LDPC codes constructed 
in Type-I are free of cycles of length four for a prime length 
n. 

Proof: Consider the block-column indexed by n — j for 
1 < .7 < tt. — 1 and let and r\ be two different block-rows 
for 1 < r.i, r' { < (8 — 1). Assume by contradiction that we have 

A(a ri ( n - j )) = A(a r> >-^). Thus n(n-j) mod n = r' t (n- 
j) mod n or n(fi — r ? ') mod n = (r» — r' t )j mod n = 0. 
This contradicts the assumption that n > j > 1 and 7^ r[. 



< n < q r ' 



1 over F g with designed distance 8 in 
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Hence primitive BCH codes of composite length n can not be 
used to derive LDPC codes that are cycles-free of length four 
using our construction. 

The proof of the following lemma is straight forward by 
exchanging, adding, and permuting a block-row. 

Lemma 8: Let (. . . , L?, . . .) be a vector of length fj, that 
has 1 at position £. Under the cyclic shift, the following two 
blocks h a and /if, of size /j x /j are equivalent, where h a 
and hi, are generated by the rows ( 1 ... U ... ) and 



1 



1, 



and their cyclic shifts, respectively. 



One might imagine that the rank of the parity check matrix 
H in ([Tol l is given by (<5— since rows of every block-row h a 
is linearly independent. A computer program has been written 
to check the exact formula and then we drove a formula to 
give the rank of the matrix H. 

Theorem 9: Let n be a prime in the range gL m / 2 J < n < 
fi = q' n — 1 and 5 be an integer in the range 2 < 8 < n for 
some prime power q and m = ord q (n). The rank of the parity 
check matrix H given by 



H 



A° A 1 A 2 
A A 2 A 4 



A° A 



5-1 



A 6 - 1 ) 



■ ■ A n ~ l 

. . _42(n-2) 

A(a l ). 



(10) 



is (<J-1>- (S-2), where A 1 

Proof: The proof of this theorem can be shown by 
mathematical induction for 1,2,..., 5 < n. We know that 
every block-row is linearly independent. 

i) Case i. Let 8 = 2, the statement is true since ever block- 
row has only 1 in every column, the first n columns 
represent the identity matrix. 

ii) Case ii-1. Assume the statement is true for 8 — 2. In this 
case, the matrix G has a full rank given by (<5 — 2)/i — 
(8 — 3). So, we have 



/ h n 


hi2 


hi 3 


hi n 


\ 





h 22 


h 23 


h 2n 










h 33 


h 3n 










: 


hin 




V o 





■•• fy(5-2)((5-2) ■ 


■ ■ fy(5-2)n 


/ 



The elements h' u s have l's in the diagonal and zeros 
everywhere using simple Gauss elimination method and 
Lemma [8] 

iii) Case iii-1. We can form the sub-matrix H2 of size (8 — 
1)/jx (8 — by adding one block-row to the matrix G. 
The last block-row is generated by 

(A(a°), ^(a 5 - 1 ), A(a 2(5 - 1} ), . . . , A(a n - 1{S - 1) )). 

All /.t — 1 rows of the last block-row are linearly indepen- 
dent and can not be generated from the previous 8 — 2 
blocks-row. Now, in order to obtain the last row-block 
to be zero at positions h {s _ 1)1 , h {5 _ 1)2 , h(5-i)(6-2), 
we can add the element hjj to the element h/g_i\j. In 



addition, the last row (row indexed by (8 — of block- 
row (5—1 can be generated by adding all elements of the 
first block-row to the first fi— 1 rows of the last block-row. 



/ hn 


hi2 


hi 3 


hin 


\ 





h 2 2 


h 23 


h 2n 










h 33 


h 3n 










: 


hin 




V 







■■ fy<5-l)n 


/ 



Therefore, the matrix G has rank of (6 — 2)/U — (5 — 3) + 
fj, — 1 = (5 — l)fi — (6 — 2). We notice that the matrix 
H has the same rank as the matrix G, hence the proof is 
completed. 

■ 

The proof can also be shown by dropping the last row of every 
block-row except at the last row in the first block-row. Hence, 
the remaining matrix has a full rank. 

Obtaining a formula for rank of the parity check matrix H 
allows us to compute rate of the constructed LDPC codes. 
Now, we can deduce the relationship between nonprimi- 
tive narrow-sense BCH codes and LDPC codes constructed 
in Type-I. 

Theorem 10 (LDPC-BCH Theorem): Let n be a prime and 
q be a power of a prime, such that gcd(n, q) = 1 and grL m / 2 J < 
77 < q m — 1, where 777 = ord n (q). A nonprimitive narrow- 
sense BCH code with parameters [n, k, d m i n ]q gives a (S — 
1, 77) LDPC code with rate (n/Lt - [(6 — l)fx — (8 - 2)])/n/i, 
where k = n — m\(S — 1)(1 — l/q)~\ and 2 < 8 < S ma x- The 
constructed codes are free of cycles with length four. 

Proof: By Type-I construction of LDPC codes derived 
from nonprimitive BCH codes using Equation ( fTOb , we know 
that every element a 1 in H^ch is a circulant matrix A(a l ) 
in H. Therefore, there is a parity check matrix H with size 
(8 — x n/Lt. H has a row weight of 77 and a column weight 
of 8 — 1. Hence, the null space of the matrix H defines an 
LDPC code with the given rate using Lemma [9] 

The constructed code is free of cycles of length four, 
because the matrix Hbch has no two rows with the same value 
in the same column, except in the first column. Hence, the 
matrix H has, at most, one position in common between two 
rows due to circulant property and Lemma [7] Consequently, 
they have a Tanner graph with girth greater than or equal to 
six. ■ 

Based on Type-I construction of regular LDPC codes, we 
notice that every variable node has a degree <5 — 1 and every 
check nodes has a degree n. Also, the maximum number of 
columns that do not have one in common is 77. Therefore, the 
following Lemma counts the stopping distance of the Tanner 
graph defined by H. 

Lemma 11: The cardinality of the smallest stopping set of 
the Tanner graph of Type-I construction of regular LDPC 
codes is fi + 1. 

Example 12: Let 77, = /i = q' m — 1, with m = 7 and q = 2. 
Consider a BCH code with 8 = 5 and length 77. Assume a to 
be a primitive element in F q m . The matrix H can be written 



5 



as 



H 



( 


1 


a 


a 2 . 


. a 126 




1 


n 2 


a 4 . 


. a 125 




1 


a :i 


a 6 . 


. a 124 


\ 


1 


a 4 


a 8 . 


. a 123 



(11) 



and the matrix H has size 508 x 16129. Therefore, we 
constructed a (4, 127) regular LDPC with a rate of 123/127, 
see Fig. Q] 

TABLE I 

Parameters of LDPC codes derived from NP BCH codes 



q 




BCH Codes 


LDPC code 
size of H 


rank of H 


2 


31 


[23, 12,4] 


(93,713) 


91 


3 


26 


[23, 12, 5] 


(104,598) 


101 


2 


31 


[31,26,3] 


(62,961) 


61 


2 


31 


[31,21,5] 


(124,961,) 


121 


2 


31 


[31,26,6] 


(155, 961) 


151 


2 


31 


[31,16,7] 


(186,961) 


181 


2 


63 


[47, 24, 4] 


(189 ,1961) 


187 


2 


63 


[61,21,6] 


(315, 3843) 


311 


2 


63 


[61,11,10] 


(567,3843) 


559 


2 


127 


[127, 113, 15] 


(1778,16129) 


1765 


2 


127 


[127, 103, 25] 


(3048,16129) 


3025 



IV. LDPC Codes Based on Cyclotomic Cosets 

In this section we will construct regular LDPC codes based 
on the structure of cyclotomic cosets. Assume that we use the 
same notation as shown in Section [TT] Let C x be a cyclotomic 
coset modulo prime integer n, defined as C x = {xq l mod n \ 
i G Z, 1 < x < n}. We can also define the location vector y 
of a cyclotomic coset C x , instead of the location vector z of 
an element a 1 . 

Definition 13: The location vector y(C x ) defined over a 
cyclotomic coset C x is the vector y(C x ) = (zq, zi, . . . , z n ), 
where all positions are zeros except at positions corresponding 
to elements of C x . 

Let I be the number of different cyclotomic cosets C x 's 
that are used to construct the matrices Hq.'s. We can index 
the i location vectors corresponding to C Xl , C X2 , . . . , C Xe , as 
y 4 ,y 2 , . . . , y e . Let y 1 { r )C x ) be the cyclic shift of y 1 (C K ) where 
every element in C x is incremented by 1 . 



A. Type-II Construction 

We construct the matrix from the cyclotomic C x as 



( y\c x ) \ 
f{ic x ) 



\ y l {i n - x C x ) j 



where y 1 ('~f : ' +1 C x ) is the cyclic shift of y l {^C x ) for < j < 
n — 1. 

From Lemma [5] we know that all cyclotomic cosets C x s 
have a size of m if 1 < x < nq^ m / 2 ^ / (q m — 1). 



We can generate all rows of Hc x , by shifting the first row 
one position to the right. Our construction of the matrix H l c 
has the following restrictions. 

• Let x < Q(y/n), this will guarantee that all cyclotomic 
cosets have the same size m. 

• Any two rows of H l c have only one nonzero position in 
common. 

• Every row (column) in H l c has a weight of m. 

We can construct the matrix H from different cyclotomic 
cosets as follows. 



H 



[ Hb 
( 



Hc 3 

y x (Ci) 



y 2 (c 2 ) 
y 2 (7C 2 ) 



(13) 



\ y 1 (7- 1 Ci) 



r(r-- 1 c 2 ) ... y%" 

where we choose the number I of different sub-matrices Hcj ■ 
The n x (i * n) matrix H constructed in Type-II has the 
following properties. 

i) Every column has a weight of m and every row has a 
weight of m*i, where I is the number of matrices H' c s. 

ii) For a large n, the matrix H is a sparse low-density parity 
check matrix. 

We can also show that the null space of the matrix H defines an 
(m,m£) LDPC code with rate ((. — l)/£. Clearly, an increase 
in t, increases the rate of the code. 

Since all cyclotomic cosets C Xl , C X2 , . . . , C xe used to con- 
struct H are different, then the first column in each sub-matrix 
H J C is different from the first column in all sub-matrices H l c 
for j 7^ i and 1 < i < I. Now, we can give a lower bound in 
the stopping distance of Type-II LDPC codes. 

Lemma 14: The stopping distance of LDPC codes, that are 
in Type-II construction, is at least I + 1. 

One can improve this bound, by counting the number of 
columns in each sub-matrix H % c that do not have one in 
common in addition to all columns in the other sub-matrices. 

Example 15: Consider n = q m — 1 with m = 5, q = 2, 
and 6 — 5. We can compute the cyclotomic cosets C\, C3 
and C 5 as d = {1,2,4,8,16}, C 3 = {3,6,12,24,17} and 
C 5 = {5, 10, 20, 9, 18}. The matrices H^, #c 3 and H c & can 
be defined based on C\, C3 and Cg, respectively. 

/1101 0001 0000 0001 0000 0000 0000 000\ 

0110 1000 1000 0000 1000 0000 0000 000 

0011 0100 0100 0000 0100 0000 0000 000 

0001 1010 0010 0000 0010 0000 0000 000 

0000 1101 0001 0000 0001 0000 0000 000 



0100 0100 0000 0100 0000 0000 0000 011 
V1010 0010 0000 0010 0000 0000 0000 001/ 



(14) 



(12) The matrix H of size (31,93) is given by 



H 



H 3 



(15) 



therefore, the null space of H defines an (5,15) LDPC code 
with parameters (62, 93). 

We note that Type-I and Type-II constructions can be used 
to derive quantum codes, if the parity check matrix H is 
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- 4,31 )-LDPC from BCH code with size (124,961) 

- (4,23)-LDPC from BCH code with size (124,713) 

- (5,61 J-LDPC from BCH code with size (315,3843) 



Fig. 1. Type I: The error performance of an (4,31) LDPC code with rate 
27/31 and H matrix with size (124, 961) based on a BCH code. 



modified to be self-orthogonal or using the nested propery of 
LDPC-BCH codes. Recall that quantum error-correcting codes 
over F q can be constructed from self-orthogonal classical 
codes over F g and F q 2, see for example [2], [3], [5], [10], 
[15] and references therein. 

V. Simulation Results 

We simulated the performance of the constructed codes 
using standard iterative decoding algorithms. Fig. Q] shows the 
BER curve for an (4,31) LDPC code Type I with a length of 
961, dimension of 837, and number of iterations of 50. This 
performance can also be improved for various lengths and the 
designed distance of BCH codes. The performance of these 
constructed codes can be improved for large code length in 
comparison to other LDPC codes constructed in [12], [13]. 
As shown in Fig. Q] at the 10~ 4 BER, the code performs at 
5.5 Eb/No(dB), which is 1.7 units from the Shannon limit. 

VI. Conclusion 

We introduced two families of regular LDPC codes based 
on nonprimitive narrow-sense BCH codes and structures of 
cyclotomic cosets. We gave a systematic method to write every 
element in the parity check matrix of BCH codes as vector 
of length [i. We demonstrated that these constructed codes 
have high rates and a uniform structure that made it easy 
to compute their dimensions, stopping distance, and bound 
their minimum distance. Furthermore, one can use standard 
iterative decoding algorithms to decode these codes. One can 
easily derive irregular LDPC codes based on these codes and 
possibly increase performance of the iterative decoding. Also, 
in future research, these constructed codes can be used to 
derive quantum LDPC error-correcting codes. 
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"Accurate reckoning: The entrance into knowledge of all 
existing things and all obscure secrets." Foundation of true 
science, Ahmes, Anc. EG. Scribe, 2000 BC. S.A.A. confirms 
that the simple work accomplished in this paper is based on 
accurate counting. 
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